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1 .0 INTRODUCTION 


This first monthly progress report details the mathematical develop- 
ments necessary for the extended. Kalman filtering and modified Bryson- 
Frazier smoothing algorithms. This is an expansion on the material 
presented at the initial meeting on 3 May 1983 including the current 
mathematical developments . 

In this initial monthly progress report, only the first twelve system 
state variables are presented with the necessary mathematical developments 
for incorporating them into the filter /smoother algorithm. Other state 
variables, i.e., aerodynamic coefficients can be easily incorporated into 
the estimation algorithm, representing uncertain parameters, but for 
initial checkout purposes are treated as known quantities. 

An approach for incorporating the NASA propulsion predictive model 
results into the optimal estimation algorithm has been identified. This 
approach utilizes numerical derivatives and nominal predictions within 
the algorithm with global iterations of the algorithm . The iterative 
process is terminated when the quality of the estimates provided no longer 
significantly improves. 


2.0 FILTERING AND SMOOTHING ALGORITHt-l 


The Space Shuttle Parameter Estimation Program utilizes optimal 
estimation techniques to provide estimates of the propulsion system 
parameters. The technique selected is the extended Kalman filter and 
the modified Bryson -Frazier smoother. By modeling the propulsion system 
parameters as time correlated random variafolos , improved estimates of 
these parameters are obtained and ai'e properly time phased by removing 
the filter induced lag by using the combined filter /smoother. The 
smoother also provides improved estimates of the initial state ;:stimates 
The system, in state-space notation, is modeled as the continuous 
dynamical system disturbed by additive Gaussian \diiite noise 

X = ^(x(t), t) + G(t) w(t) -’r u(t), x(o) = ^ (1) 

where 

X = n-dimensional state vector 

X = Gaussian initial condition vector with covariance P 
— o o 

w(t) = p-dimensional white, zero-mean white Gaussian noise with 
covariance 

E[w(t) w'^(t)] = Q(t) 6(t - t) 
u(t) = n-dimensional control vector. 

The elements of the vector x(t) represent vehicle position, velocities 
attitudes, angular rates, aerodynamic and propulsion parameters, measure 
ment biases, etc- Elements of u_(t) include known control inputs such 
as SSME power level commands . 



The system described by equation (1) is observed at discrete times, t^ 
with not all states being directly measured. Some measurements are non- 
linear functions of the elements of the state vector x(t) . In general 
the measurement process is described as 

2^ = h^(x(t|^)) -f 4 (2) 

where 


Zj^ = m-dimensional observation vector 

k = functional representation of the measurements in terms of 


the states 


Vj^ = m-dimensional, zero-mean, while Gaussian noise sequence with 
covariance 


eEv. 


-a 


R. 

1 


6 . . 

3-0 


Examples of the elements of the observation vector z^^ include radar 
measurements of range, azimuth, and elevation from the radar site to the 
vehicle . 

It is assiimed that the system process noise vector w(t) and the 
measurement noise vector v are uncorrelated. Also, the system state 

"tC 

initial condition vector x is not correlated with either of these two 

—o 

noise vectors . Therefore 


E[w(t) v^] = 0, E[w(t) x"^] = 0, E[x v^] = 0 

“ “K ‘~0 —O “K 

wh'iSre the superscript T denotes transpose. For later reference, the 


follovving matricies are defined 


cf_(x(t), t) 
FMt>. t) = ■— ijTtj — 


( 3 ) 


and 

0h(x(t ) ) 

These matricies are linearizations of the dynamics and measurement models 
respectively I evaluated eibout either a nominal or reference value of the 
state, or about the sttite estimate. 

2 .1 E >ctended Kalman Filter Algorithm 

The extended Kalman filter algorithm is in essence a conventional 
linear Kalman filter algorithm applied to a mathematical model resulting 
from the linearization of the system model equation (1), and measurement 
process, equation (2), about a current state estimate. The filter yields 
optimal estimates if the linearization is accurate, i.e., the state esti- 
mate closely approximates the true state. The derivation cif the algorithm 
can be found in reference [13. 

The algorithm proceeds as follows. After initialization of the state 
estimate and covariance, the state estimate and covariance are propagated 
forward in time until a measurement update is available, by 

X = ^(x(t) , t) , \ 1 ^ ~ \ ^ 

and 

P(t) = P(x(t), t) P(t) + P(t) F(x(t), t)^ + G(t) Q(t) G(t)^ 


( 6 ) 


At the measurement time, the state estir^ate and covariance are updated 
by 

4 <+) = 4 '-) 

and 

P^( + ) = (I - Kj^ 

where the (-) and (+) represent the appropriate values just before 
and just after the update. The updated values are used to reinitialize 
the time propagation equations (3) and (4) for integrating up to the next 
measuremerit time . The Kalman gain matrix is computed a,s 

This algorithm is repeated until the last time point, t^^, is processed. 

For later use in the smoother algorithm, various combinations of the state 
estimates (x), measurements (z), linearized dynamics matrix (F) and 
measurement matrix (H), measurement noise covariance (R) and estimation 
error covariance matrix (P) must be stored for each time instant to be 
processed by the smoother algorithm. 

2 .2 Modified Bryson -Frazier Onoother Algorithm 

Ibe operation of the smoother algorithm is similar to the filter 
algoritfii". except in reverse time. The derivation of this smoother algorithm 
• is' found in reference [2]. This fixed interval smoothing algorithm pro- 
vides optimal estimates given all the measurements in comparison to the 
filtering algorithm providing optimal estimates given the previous 
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measurements processed. Therefore the smoother provides improved estimates 
in addition to removing the time lag induced by the filter algorithm. 

The smoothing algorithm adjoint variables, ^ and A a’ a "initialized" 
at the final time processed by the filter, T, 


UT-) = 6,^ 


N,T 


and 


( 10 ) 


A(T-) = H. + R J"'' H , 6^ 

N N N N N N t 


N,T 


(11 ) 


If T is not an observation time, ^ and A are zero. The adjoint 
variables are propagated in reverse time to the next previous measurement 
time by 


A = - F(x(t), t) X , \ (12) 

A = - F(x(t), t)*^' A - A F(x(t), t) i\^) 

At the time of an available measurement, t , the adjoint variables 

K> 

are updated by 

X(_) = x(+) - H^(H^ \ ^ 

and 

. A(-) = (I - A( + ) (I - I^) + (15) 
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The smoother state estimate and error covariance fitre obtained using the 
filter estimate and covariance and the adjoint variables by 

X (t) c “ P(t) y^(t) (16) 

and 

p'^’Ct) « P(t) - P(t) A(t) P(t)» (17) 


Due to the potential number of time points to be processed, smoother 
estimates may only be computed at the discrete measurement times . For 
this approach the propagation equations (10) and (11) are replaced by 




(18) 


and 


A (+) = O A (-) $ 
k k+1 ^ ^ k 


(19) 


where is the state transition matrix formed with the linearized 

k 

dynamics matrix F to propagate the adjoint variable from time 
to time t|^. The algorithm continues in reverse time until the initial 
time is reached. 


2 .3 Iterations with the Filter/Snoother Algorithm 

The performance of the filter /smoother algorithm is a direct result 
of the accuracy of \,aq linearization. Repeated operations of the algorithms 
•with adjustments in initial state estimates and covariance in each cycle 
can yield improved estimates.. This technique is known as global iterated 



filtering as defined in reference [3]. Each cycle of operating the 
algorithms would yield increasing imps-t' vements in the state estimates . 

This feature of the algorithm operation is of special interest to 
the propulsion parameter estimation problem using the NASA predictive 
models. Initial, or nominal, values of the pax*ameters of interest can 
be used to obtain the necessary partial derivatives indicated earlier „ 
From operating the algoritiim improved estimates of those pai*ameters are 
obtained. Using these improved estimates, more accurate partial deriva- 
tives are obtained for use in the algorithms. This process is continued 
until there is in essence no change in the partial derivatives or quality 
of the state estimates. If the linearization is accurate, the measure- 
ment residual should be a white noise process with known covariance . 
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3.0 FILTER/SMQOTHSR ALGORITHM SYSTEM AND MEASUREMENT MODEL 

Hie usefulness of the filter/smoother algorithm is to provide esti- 
mates of the system states from the observed motion and dynamics while 
the system is driven by known and unknown elements. These unknown elements 
are elements of the system state vector to be estimated. The evolution of 
motion resulting from these known and unknown elements is assumed to be 
suitably represented for this study by a six degree -of -freedom (6 DOF) rigid 
body equations of motion . These equations are presented and discussed in 
section 3.1 . 

To implement these equations i,ito the filter /smoother algorithm 
presented in section 2.0, a linearization of the system state and measure- 
ment models is required. These linearized equations are presented in 
section 3.2. 


3.1 Equations of Motion and Measurement Equations 


3.1.1 Rigid Body Equations of Motion 

The rate of change of vehicle velocity in body coordinates , v 


(B) 


a& a result of external forces acting on the vehicle is described by 
2 


.(B) 


pAv 


2m 


-(B) (B) 

m B I (I) , (I), (B) ^ ^ 


m 


m 


( 20 ) 


where 

p = atmospheric density 

A = aerodynamic coefficient referenced area 

v^ = magnitude of vehicle velocity relative to the surrounding air 


mass 




m " vehicle mass 
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» aerodynamic force coefficient vector 
(1) (I) 

2 _ (r ) * gravity vector in inertial coordinates 

B angular rotation of the body relative to the inertial frame 
{ B ) 

B resultant thrust force vector in body coordinates 
CB ) 

f „ » resultant plume force vector in body coordinates 

(I ) 

The rate of change of vehicle position in inertial coordinates , r > 
is then obtained by 

^(I) ^ ^(B) (2- 

I B 

where C is the transformation matrix from body coordinates to 

I B 

inertial coordinates . The elements of the C transformation matrix 
are obtained from the resulting Euler angles defined by 


<p 


1 sin 9 tan 0 cos(ptaii0 


P 

• 

e 


0 cos(p -sincp 


q 



0 sin(psec0 costpsec0 


r 


where <p, 0, and are roll, pitch and yaw attitudes respectively. 

The roll , pitch and yaw rates of the body relative to inertial coordinates 
are p, q, and r respectively. Finally, the rate of change of the body 
rates relative to inertial is given by 
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to 


. PAv 

- Cl] + 

2 — m 


2 


) X 

—eg -£ 


- (0 X (Tto) + 

— — — T p 


(23) 


where 


1 a vehicle moments of inertia matrix 

= aerodynamic coefficient referenced length and moment 

coefficient vector 
(b) 

r, a vehicle center-of -gravity vector in body coordinates 

“v-g 



■h aerodynamic coefficient reference position in body coordinates 


= resultant thrust torque vector in body coordinates 
( B ) 

Tp = resultant plane torque vector in body coordinates 


The equations of motion represent the first twelve elements of the system 
state vector. These equations are summarized in Table 3.1 .1-1. 

The moment of inertia matrix I in general is given by 



-I 

-I 

X 

xy 

zx 


I 

-I 

xy 

y 

yz 


-I 

I 

zx 

yi. 

z 


(24) 


for the moment axis terms* i.e., 

X . e . ^ X . 

zx 



and the product of inertia terms, 


TABLE 3.1 
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The aerodynamic force and moment coefficients and plume forces are 
defined as functions of angle-of -attack, a, and angle-of -sideslip, 3, as 
shovn in Figure 3.1 .1-1. The body referenced relative velocity vector, 
removing the wind velocity, v , from the vehicle velocity, is given by 


B) (B) B LL (LL 

v=v - Cv=v - G V 
—r— -w — -w 


(25 ) 


where v 
— w 


(LL) 


is the local -level referenced wind velocity vector. The 


following equations define ct and 3 in terms of the components of v 

*~ir 


V 

-1 ^3 

c =5 tan (— — ) 


(26) 


-1 2 

3 = sin (- — ) 

V 

m 


(27) 


where 


, 2 2 2 , 

V = (v + V + V ) 

m r^ r^ r3 


(28) 


The resultant thrust force f„ 

-T 


(B) 


- X 

-T 1=1 1 


T. 

X 

0 

0 


A 

= .2 
1=1 


is expanded as 


B ^ ) 

C f ^ 

i 


(29) 


B 

where the transformation matrix transforms the magnitude of thrust 

for each thrusting device, f^ , from its center -line to the body 

i 

coordinates . The general eqaution for f ^ is 

i 


UI >. 

SE oc 

HOO 

<OI- 

J-IU 

UJUJ Ul 

tc>> 



FIGURE 3.1 .1-1 - Body Referenced Wind Axis System 
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= f 


T. 

ivac 



A 


e 


where 


= vacuum thrust 

ivac 

P = motor exit pressure 
i 

A = motor exit cone area 
e 

B ^ 

The matrix is different for the SSME’s and SRB's and is given by 



B MP MP G, G 


B 4 
C 


SSME 


SRB 


( 30 ) 


Where 

B MP 

C = transformation from the engine mount plane to the body 
coordinates 

MP G 

C = transformation from the gimbal reference plane to mount plane 
(structural deformation) 

G 

C = txansforraation from enterline to the gimbal reference plane 
B ^ 

C = transformation from SRB nozzle centerline to the body 
coordinates (gimbal angles). 
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The resultant thrust torque is the summati.on of the torque contribution 
from each thrusting device and is given by 




where 


2 (r<®> 

ill 


B), 

r ) C. 
-eg 1 


■T. 

.1 

0 

0 


(31 ) 




(B) 

I 

i 


= body coordinates of the thrust reference point for the 
i' thrusting device. 


3.1.2 Measurement Equations 

The measurements assumed available for the filter /smoother algorithm 
include inertial platform acceleration and attitudes , ground based radar 
tracking, SRB' s head pressure, SSME' s chamber pressures, liquid flow 
rates, pressurant flow rates. The ET volumetric levels are available; 
however, due to their limited number (4), they may only be used for 
alternate checks of the filter /smoother algorithm performance. 

The propulsion related measuremetns will be treated in a separate 
section. In the following, the inertial platform acceleration measurements, 
attitude measurements and ground based tracking measurements models will 
be described for later linearization. 
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3. 1.2.1 Platform acceleration Measurements 

Accelerometers mounted orthogonally on an inertially stabilized 
platform, not located at the vehicle center of gravity, sense externally 
applied special forces and accelerations due to body rotation. Ihe 
accelerometer measurement is modeled by 


“tn 


S.,P P^P' 
= C C 


pAy 


m 


2m 




j(B) (B) 

ijT -P 
+ 


m 


m 


+ 0) X 0) X 




-r«») 

~cg 


+ U) X 




-eg 




(s) 

+ Zo 
—a 


(32 


where 

S P 

C = transformation from platform coordinates to sensing coordinates 
P P' 

C = transformation from misaligned platform coordinates to 
platform coordinates 

P * B 

C = transformation from body to misaligned platform coordinates 

(b ) 

r ' = body coordinates of the platform center 

(s ) 

b = accelerometer bias vector 
—a 

v^^^ = accelerometer measurement noise vector 
—a 

3. 1.2. 2 Placform Attitude Measurements 

The inertially stabilized platform for the STS is a four axis IMU 
with a redundant roll axis [4] . Vehicle body attitudes are generated 
via guaternions [5]. It is assumed that an equivalent representation 
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can be made to obtain vehicle attitude by a three rotation sequence of 
roll, pitch, yaw to transfoi'm from inertial to body coordiantes . This 
approach has been used in reference [61 .. 

The attitude angle measurement model is given by 


g(S) 


= 0 


4 . 

+ ^ 




(S) 


(33) 


where 


P P* 

^ = platform misalignment bias vector (used to formulate C ) 

( S ) 

Vq =: attitude measurement noise vector. 

The transformation mafa’ix used to transform from body to inertial 
coordinates in terras of the elements of the 0 vector is given by 


cos0cos<|< sin<psin0cos(|j 

-costpsineji 


cos<psin0cos(j) 

+sin<psint(> 


I B 
C 


cos0sin(|> 


sin(psin0sin<)> 

+COS(pCOS(|j 


cos(psin0sin(^ 

-sin(pcos(|j 


“Sin0 


sincpcos0 


CQS(pcos0 


(34) 


3. 1.2. 3 Ground Based Tracking Measurements 

Ground based radar tracking devices can provide measurements of 
range, azimuth and elevation from the radar sight to the vehicle. Azimuth 
and elevation are established relative to the sight's local level. If 


ORJG^NAt 

OF POOR 

the tracking device is a passive optical tracker (not laser) then only 
azimuth and elevation measurements are available requiring more than one 
to establish position information. 

Defining Xj y, and z as the local east, north and up position of the 
vehicle relative to the ground based tracking device, the radar measurement 
equations are given by 


lei 

qUAUVTV 


p = (x^ + y^ + z^) + b + V 

P P 

A = tan“^ (— ) 4- b + v 
y A A 


( 35 ) 

( 36 ) 


E s= tan (- 




2 2 
X + y 


) b + A + V 
E E E 


where 


( 37 ) 


b , b_, b„ 5= range, azimuth, elevation biases 

p A hi 

AE = atmospheric retraction correction 


, V , V , v_ = range, azimuth, elevation measurement noise. 

p A hi 


The position vector of the vehicle relative to the tracking device is 
given by 


X 

A . (LL) LL ECF .ECF ECI (I) (ECF) 
y =Ar - C [ C r -r 


( 38 ) 
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where 

LLcECP ^ transformation from earth center fixed to local level 

ECF^ECI _ Qa]-th centen)\d inertial to earth centered fixed 
(ECF ) 

“ position vector of tracking device in ECF coordiantes . 

““vDR 

The transformation matrix is given by 


-sink 

“SinLcosX 

cosLcosX 

cosX 

-sinLsinX 

cosLsinX 

0 

cosL 

sinL 1 


J 


( 39 ) 


where L and \ are the geodetic latitude and east longitude of the 

ECF ECl 

device. The transformation matrix C is given by 


ECF^ECI 

C 


cosCm^Ct - sin[m^(t 

-sin[m^(t - t^^p)] cos[(Ug(t 
0 0 






ERNP] 


( 40 ) 


Where 


0 ) = earth rotation rate 

^RNP ~ 'tinie tag for RNP matrix 


The position vector, r' _ , of the tracking device is given by 

•^KDK 


'! 
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( - 

V cos^L + 


7 ! ^2 27 

(1 - e) sin L 


+ h) cosLcoc\ 


(ECF) 

-RDR 


( — - ' ' - - ' + h) cosLsinX 

V cos^L + (1 ” e)^ sin^L 


Rg(1 - e)‘ 




+ h) sinL 


cos^L + (1 - e)^ sin^L 


(41 ) 


where 


R E equatorial radius of Fisher ellipsoid 
E 

e = flattening of Fisher ellipsoid 
h = altitude of the device above Fisher ellipsoid 



3 .2 Linearized System State and Measurement Equations 


The vehicle equations of motion are nonlinear functions of their 
motion variables and are implicit functions of other elements of the 
system states . The measurement equations involve similar* function rela- 
tionships . The linearizations for the filter /smoother algorithm require 
partial derivatives with respect to the motion variables, i.e., v ' and 
0 , and with respect to other elements of cha state vector . yielding 
explicit functional relationships for the elements of interest. 

For system state equations the partial derivatives wili be presented 
in section 3.2.1 for the state elements in order of occurrenctj. for the 
first twelve states. Other partial derivatives for candidate state ele- 
ments will follow in section 3.3.1. The measurement equation partial 
derivatives for the first twelve states will be presented in section 3.2.2. 
Partial derivatives of the measurement equations for other candidate states 
will be presented in section 3 .3 .2 . 

The resulting partial derivatives are imbedded into the linearized 
system state matrix, F(x(t), t), as shown in Figure 3.2-1 . A corresponding 

linearized measurement matrix, H(x^), is similarly formed with the 

i\ 

measurement equations' partial derivatives. 

3.2.1 System State Partial Derivatives 

Partial derivatives of each of the equations listed in Table 
are developed in their order of occurrence with respect to the order of 
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tho corresponding rtates. Partial derivatives of thrust terms are presented 
as though for a single device. 
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Inertial Position Rate Equation 

* ( X ) 

Ihe first nonzero partial derivative of the r equation is with 
(B) 

respect to v : 


3 

a^(B) 




V. 


(42) 


The second nonzero partial derivative is with respect to This partial 

derivative results in a third order tensor and occurs frequently in later 
developments. The generalized form is presented in Appendix A. 


Body Velocity Rate Equation 

The partial derivative of v' with respect to r' ' is given for 


altitude terms approximately as 

.T 


3r 


(I) 


a_ 

6h 


(I) 


.(I) 


(43) 


where 


0V 


(B) 


3h 


Av - pAv 

m 3p ^ ' 

3m 3h ^ 


3v pAvf ac^ pAvf 3c- 

m m_ -€ oa ^ rn_ -f op 

m ah 2m 3a 3h 2m SP 3h 


(B) 


+ -( c 

m 3p 


dp 3f 
s -p 

3h 3h 


3f^®^ 

-P 


s 


3g 
3h 3h 




3f^°^ Aft 

33 3h- 


( 44 ) 
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The partial derivatives of 


• 8a 86 Sa 86, 

8v'®> ’ 8v<®> ' 8v<®> ’ 


8v 


and 


m 


3h 


occur frequently and are given in Appendix B. 

Dl (l) (l) 

The gravity vector C ) partial derivative with respect 

4 . (I) . 

to r xs 


B I 
C 


ft (I)/ (I)\ 

dg (r ) 


3r 


(I) 




(I) |3 


,^1^2 


,^1^3 


- 1 




,^ 2^3 


- 1 




,^2^3 


where 


M = gravitational constant. 

(B) 


3v 


The partial derivative, 


3r 


(I; 


, is the sum of the matxicies in equations 


and 


•(B) (B) 

The partial derivative of v with respect to v is given by 


3v 


(B) 


3v 


(B) 


pAv 3v 
m m 

— 
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» skew symmetric matrix made from the elements of the vector w 
and equivalent to the cross product operator w x ( ) , 

• (B) , 

The partial derivative of v with respect to 0^ is 
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The 


peirtial derivatives of are given in Appendix B and the partial 


3v 

~r 


derivative is given in Appendix A. The last partial derivative 

OU 


is given in Appendix C. 
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The partial derivative of v with respect to « is 


3v 


(B) 


pAv 


m 


3u) 


CB), 

ar-+ ^ 


2m 


(48) 


original PAGlS m 

OF POOR QUALITY 


Euler Angle Rate Equation 

The Euler angle rate equation is a function of both the Euler angles 
and the inertial rates. The linearization will yield the two associated 
matricies . 

First with respect to the vector 0, the following matrix results 


39 
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The partial derivative of 


0 with respect to to is 
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Inertial Angular Acceleration Equation 

The first partial derivative of this equation is with respect to the 
vector Using the approximation indicated in equation , this 


partial derivative is 
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Next, with respect to the vector v , the partial derivative is 
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The partial derivative with respect to the vector 0 is 


= tir' ^ X pv^A 


2 0 
pv M 0c , Dv A 0c 

+ (— =! 21 + - r<®>) X I" 1“ 

2 0a ^ 2 0a^ 00 


2 9 

pv Ad 0c /r^v /. ovA 0c 

+ i-2 2 L+ (r<B) . ^(B), P'm ^ °Sf as 

2 03 + '-A ^ ^ 2 03^ 00 


0T ^ 0T 
+ ~P . •"•P 9& 1 

0a 00 03 00 


The final partial derivative for the first twelve states is with respect 
to the vector w. This operation yields 
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3.2.2 Measurement Partial Derivatives 

The measurements assumed to be available, as discussed earlier, include 
ground based radar tracking, inertially stabilized platform attitudes 
relative to the vehicle body, and stabilized platform mounted 3 axis 
orthogonal accelerations. As with the state dynamics matrix, the measure- 
ment equations are linearized aobut the best state estiarates- 


Radar Track Measurement Equation 

Referring to the radar track measurement equations, the required 
partial derivatives are 
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The last partial derivative in each of these equations, ■ — - T rT ' " > is 
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The rest of the required partial derivatives are 
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Inertially Stabilized Platform Attitude Equation 

The inertial platform is assumed to provide attitude angle measure- 
ments of the true attitude plus an attitude bias plus measurment noise . 

The partial derivative of the measured attitudes with respect to the 
vector ^ yields an identity matrix. 

Accelerometer Measurement Equation 

The accelerometer senses sepcific body forces excluding gravity along 

the sensing axes. With reference to the accelerometer equation, the 

(I ) 

partial derivative with respect to r is 
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The partial derivative with respect to v yields 
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For the partial derivative of the aoeelerometer with respect to 
the vector 0 ^, the measurement equation is temporarily rewritten as 


-m -■ —a 
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where the vector s_ represents the sum of the aerodynamic, thrust, 

P' B 

plume and rotational coupling terms. The matrix C is the same 


I B 

matrix as C . The required partial derivative results from 
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The partial derivative on the right hand side is developed in Appendix A 

( B ) 

with the vector £ representing the sum of the terms indicated above . 

The final partial derivative for the accelerometer measurement is 
with respect to the body rotation vector w. Defining 
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3 .3 Propulsion Parameter States and Measurements 

A candidate approach for incorporating the NASA propulsion model's 
capabilities has been identified. This approach utilizes nominal pre- 
dicted values of thrust, pressure, propellant and pressurant mass flow 
rates, and utilizes sensitivities or partial derivatives of these 
variables with respect to the independent parameters selected for esti- 
mation by the algorithm. 

The approach is to include deviations from nominal values of 
measured chamber pressure, power level, propellant and pressurant mass 
flow rates as states . The models assumed for these deviations are time 
correlated random processes . Then as states , partial derivatives of 
the first twelve states with respect to these variables will be required. 

For the SSME and SRB, this modeling approach is discussed in the 
following. Mditionally, the necessary partial derivatives of the first 
twelve state variables with resepct to the additional states are presented. 

3 .3 .1 SSME Propulsion Parameter Model 

For the SSME, the total actual values of vacuum thrust and oxidizer 
mass flow rates are modeled by 

^T = ^T + ^^T 
nom 

and 

m^ = m + M . (69) 

2 2 2 
nom 
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The measurements of fuel mass flow ratei pressurant mass flow rates and 
power level are modeled as 


m„ e m,, Am„ + b» -f s* 

"2 "2 "2 " h , ■’ h , 
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= m 
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+ b' 
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•f s ' 


nom 


and 




( 72 ) 


PL = PL + APL + b + s . 
nom PL PL 


( 77 ) 


Ihesa measured quantities include measui-ement noise s ^ ^ and potential 
bias states b^ ^ modeled as random constants. In these measurements , 
the A'd variables are to be included as states in the estimation 
algorithm. If the nominal values are zero or unknown, then the A'd 
variables absorb the entire estimate. Where required, the estimate for 
the variables used in the estimation algorithm is formed using the 
nominal and the estimate of the deviation, etc. In example, thrust and 
fuel mass flow rate estimates are formed as 
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( 74 ) 
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Tlie deviation or A'd measurement variables are modeled as time 
correlated random variables. This permits these varaibles to vary within 
a band of frequencies. The typical model is then given as 



“ A( 
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1 



(76) 


where the parenthesis ( ) would be replaced by the variables, i.e,, m . 

* 2 
ie 

For the SSME, an additional variable Ac is modeled as in equation 
76 and included as a state mult variable with the A'd measurement 
variables . 

'The thrust deviation is expanded as in the following truncated Taylor 
series as a function of the independent parameters . 
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In the ^ and equations, with equation 77 replacing f^ , 

i 

the partial derivatives of f^ with respect to the A'd variables are 


obtained directly from equation 77. 
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It is desirable to include vehicle mass as a state. Iho SSMB's 
system contribution to the mass deviation is given by 




SSME's 


2 (Am 


H, 


+ Am^ 


- Am 


2 . 

1 


'2. 

X 


H, 



(78) 


In equation 78, the Am contribution to the mass deviation is not 

^2 

available from measurements. As with the thxuet deviation, this quantity 
is formed as 
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which is in terms of other estimated state variables . In equations 77 
and 79, the deviation in mixture ratio, AMR, is obtained algebraically 
from 


AMR = 



nom 
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3,3.2 SRB Propulsion Parflinstar Model 

The approach for the SRB modeling follows closely that used for the 
SSME, Candidate independent parameters includo propellant burn rate 
exponent, a, and motor efficiency coefficient, c^ . Others can be added 
using this technique. 

The actual value of vacuum tlirust is given by equation 68 . The 
only measurement available for the SRB is the total pressure at the 
forward head end of the motor case and is modeled as 



+ AP_ b + s 


H 


nom 




Pr 
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'H 


(81 ) 


where b^ ^ and ^ represents a bias and measurement noise respectively. 

The independent parameters, Aa and Ac , are included in the 

m 

model as states . The model assumed can be as given by equation 76 or 
another suitable dynamical process, i.e., random constant. 

The thrust deviation is given by the following truncated Taylor 
series as a function of the candidate independent parameters . 
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The partial derivatives for the v and m equations with respect to 

the independent parameters are obtained directly from equation 82 . The 

mass deviation equation for the SRB is given as 
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The head pressure deviation; AP , is expanded similarly 
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3.3.3 Vehicle Mass State Variable 

The total rate of change of vehicle mass is given by 


^^ (m) ~ '^'ssme 


"'SRB ^SSME ^SRB ''' "’nON-CONSUME 

nom nom 


The first two terms in this equation are the apriori assumed nominal 
values. The third and fourth terms were discussed earlier. Ihe last 
term should be zero. 

The equations, state and measurement, in which mass occurs are the 

♦ (B) 

V and a^ equations. Assuming equation 86 can be summarized as 

• ♦ 

m + ZAm. then the mass can be written as m + ZAm. . Replacing this 
i ^ i 

expression for the mass in the two indicated equations yields the 
following partial derivatives with respect to the individual 
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4.0 PROJECTED ACTIVITIES DURING UPCOMING MONTH 


The mathematical development of the system and measurement equations 
and their linearization will be concluded. Included will be partial 
derivatives with respect to r > v , and I for center -of -gravity > wind 

""iJj 

velocity! and moment of inertia uncertainties. Additionally j partial 
derivatives with respect to platform tilt errors in the accelerometer 
measurement equations will be formulated. 

The estimation program development will be initiated. The program 
will be structured to permit additional parameter states to be included 
with minimal impact to the program. As experience is gained with the 
operation and results of the program! changes will inevitably be required 
Initial checkouts will focus on the propulsion parameter estimation with 
other parameters assumed known, i.e., aerodynamic coefficients. Tnese 
will be included in the state model later for a more extensive evaluation 
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PARTIAL DERIVATIVE OF THE VECTOR 


I B 
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V wrt 
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This partial derivative is one of several that occurs frequently in 
the formulation of the linearized system state and measuremetn eqautions. 
The desired partial derivative is 
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The resulting matrix is given in Table A-1 . 
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APPENDIX B 
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The following equations define the last three required partial derivatives 
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APPENDIX C 

PARTIAL DERIVATIVE OP THE VECTOR V wrt 6 


Ihe third of the frequently occurring required partial derivatives 
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resulting matrix is given in Table C-1 . 
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